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Abstract. We consider the Jaynes-Cummings model of a single quantum spin s coupled to 
a harmonic oscillator in a parameter regime where the underlying classical dynamics exhibits an 
unstable equilibrium point. This state of the model is relevant to the physics of cold atom systems, 
in non-equilibrium situations obtained by fast sweeping through a Feshbach resonance. We show 
that in this integrable system with two degrees of freedom, for any initial condition close to the 
unstable point, the classical dynamics is controlled by a singularity of the focus-focus type. In 
particular, it displays the expected monodromy, which forbids the existence of global action-angle 
coordinates. Explicit calculations of the joint spectrum of conserved quantities reveal the mon- 
odromy at the quantum level, as a dislocation in the lattice of eigenvalues. We perform a detailed 
semi-classical analysis of the associated eigenstates. Whereas most of the levels are well described 
by the usual Bohr-Sommerfeld quantization rules, properly adapted to polar coordinates, we show 
how these rules are modified in the vicinity of the critical level. The spectral decomposition of the 
classically unstable state is computed, and is found to be dominated by the critical WKB states. 
This provides a useful tool to analyze the quantum dynamics starting from this particular state, 
which exhibits an aperiodic sequence of solitonic pulses with a rather well defined characteristic 
frequency. 

1 Introduction. 

The Jaynes-Cummings model was originally introduced to describe the near resonant interaction 
between a two-level atom and a quantized mode of the electromagnetic field [Ij. When the field 
is treated classically, the populations of the two levels exhibit periodic Rabi oscillations whose 
frequency is proportional to the field intensity. The full quantum treatment shows that the possible 
oscillation frequencies are quantized, and are determined by the total photon number stored in the 
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mode. For an initial coherent state of the field, residual quantum fluctuations in the photon 
number lead to a gradual blurring and a subsequent collapse of the Rabi oscillations after a finite 
time. On even longer time scales, the model predicts a revival of the oscillations, followed later 
by a second collapse, and so on. These collapses and revivals have been analyzed in detail by 
Narozhny et al [3], building on the exact solution of the Heisenberg equations of motion given by 
Ackerhalt and Rzazewski [2]. Such complex time evolution has been observed experimentally on the 
Rydberg micromaser [3]. More recently, a direct experimental evidence that the possible oscillation 
frequencies are quantized has been achieved with Rydberg atoms interacting with the small coherent 
fields stored in a cavity with a large quality factor [5]. Another interesting feature of the Jaynes- 
Cummings model is that it provides a way to prepare the field in a linear superposition of coherent 
states |6j. In the non resonant case, closely related ideas were used to measure experimentally the 
decoherence of Schrodinger cat states of the field in a cavity • 

In the present paper, we shall consider the generalization of the Jaynes-Cummings model where 
the two-level atom is replaced by a single spin s. One motivation for this is the phenomenon 
of superradiance, where a population of identical two-level atoms interacts coherently with the 
quantized electromagnetic field. As shown by Dicke [8], this phenomenon can be viewed as the 
result of a cooperative behavior, where individual atomic dipoles build up to make a macroscopic 
effective spin. In the large s limit, and for most initial conditions, a semi-classical approach is quite 
reliable. However, the corresponding classical Hamiltonian system with two degrees of freedom is 
known to exhibit an unstable equilibrium point, for a large region in its parameter space. As shown 
by Bonifacio and Preparata [9], the subsequent evolution of the system, starting from such a state, is 
dominated by quantum fluctuations as it would be for a quantum pendulum initially prepared with 
the highest possible potential energy. These authors have found that, at short times, the evolution 
of the system is almost periodic, with solitonic pulses of photons separating quieter time intervals 
where most of the energy is stored in the macroscopic spin. Because the stationary states in the 
quantum system have eigenenergies that are not strictly equidistant, this quasiperiodic behavior 
gives way, at longer times, to a rather complicated pattern, that is reminiscent of the collapses and 
revivals in the 5" = 1/2 Jaynes-Cummings model. 

An additional motivation for studying the large s Jaynes-Cummings model comes from recent 
developments in cold atom physics. It has been shown that the sign and the strength of the two 
body interaction can be tuned at will in the vicinity of a Feshbach resonance, and this has enabled 
various groups to explore the whole crossover from the Bose Einstein condensation of tightly bound 
molecules [10] to the BCS condensate of weakly bound atomic pairs [HI [12]. In the case of a fast 
sweeping of the external magnetic field through a Feshbach resonance, some coherent macroscopic 
oscillations in the molecular condensate population have been predicted theoretically [13], from a 
description of the low energy dynamics in terms of a collection of N spins 1/2 coupled to a single 
harmonic oscillator. This model has been shown to be integrable [14] by Yuzbashyan et al. who 
emphasized its connection with the original integrable Gaudin model [15J. It turns out that in the 
cross-over region, the free Fermi sea is unstable towards the formation of a pair condensate, and 
that this instability is manifested by the appearance of two pairs of conjugated complex frequencies 
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in a linear analysis. This shows that for any value of (which is also half of the total number 
of atoms), these coherent oscillations are well captured by an effective model with two degrees of 
freedom, one spin, and one oscillator. Therefore, there is a close connection between the quantum 
dynamics in the neighborhood of the classical unstable point of the Jaynes-Cummings model, and 
the evolution of a cold Fermi gas after an attractive interaction has been switched on suddenly. 

The problem of quantizing a classical system in the vicinity of an unstable equilibrium point has 
been a subject of recent interest, specially in the mathematical community. The Bohr-Sommerfeld 
quantization principle suggests that the density of states exhibits, in the h going to zero limit, 
some singularity at the energy of the classical equilibrium point, in close analogy to the Van Hove 
singularities for the energy spectrum of a quantum particle moving in a periodic potential. This 
intuitive expectation has been confirmed by rigorous [TBI flTl [18] and numerical [19] studies. In 
particular, for a critical level of a system with one degree of freedom, there are typically | In h\ 
eigenvalues in an energy interval of width proportional to h around the critical value [17^ I18j . A 
phase-space analysis of the corresponding wave-functions shows that they are concentrated along 
the classical unstable orbits which leave the critical point in the remote past and return to it in 
the remote future . To find the eigenstates requires an extension of the usual Bohr-Sommerfeld 
rules, because matching the components of the wave-function which propagate towards the critical 
point or away from it is a special fully quantum problem, which can be solved by reduction to a 
normal form |16[ [20] . 

Finally, we will show that the spin s Jaynes-Cummings model is an example of an integrable 
system for which it is impossible to define global action-angle coordinates. By the Arnold-Liouville 
theorem, classical integrability implies that phase space is foliated by n-dimensional invariant tori. 
Angle coordinates are introduced by constructing n-independent periodic Hamiltonian flows on 
these tori. Hence each invariant torus is equipped with a lattice of symplectic translations which 
act as the identity on this torus, the lattice of periods of these flows, which is equivalent to the 
data of n independent cycles on the torus. To get the angles we still must choose an origin on 
these cycles. All this can be done in a continuous way for close enough nearby tori showing the 
existence of local action angle variables. Globally, a first obstruction can come from the impossi- 
bility of choosing an origin on each torus in a consistent way. In the case of the Jaynes-Cummings 
model, this obstruction is absent (see ^1], page 702). The second obstruction to the existence of 
global action-angle variables comes from the impossibility of choosing a basis of cycles on the tori 
in a uniform way. More precisely, along each curve in the manifold of regular invariant tori, the 
lattice associated to each torus can be followed by continuity. This adiabatic process, when carried 
along a closed loop, induces an automorphism on the lattice attached to initial (and final) torus, 
which is called the monodromy [21j. Several simple dynamical systems, including the spherical 
pendulum [21] or the the champagne bottle potential [22j have been shown to exhibit such phe- 
nomenon. After quantization, classical monodromy induces topological defects such as dislocations 
in the lattice of common eigenvalues of the mutually commuting conserved operators [23]. Inter- 
esting applications have been found, specially in molecular physics [21]. In the Jaynes-Cummings 
model, the monodromy is directly associated to the unstable critical point, because it belongs to 
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a singular invariant manifold, namely a pinched torus. This implies that the set of regular tori 
is not simply connected, which allows for a non-trivial monodromy, of the so-called focus-focus 
type [25l [26l [27] . The quantization of a generic system which such singularity has been studied in 
detail by Vii Ngoc |28] . Here, we present an explicit semi-classical analysis of the common energy 
and angular momentum eigenstates in the vicinity of their critical values, which illustrates all the 
concepts just mentioned. 

In section [2] we introduce the classical Jaynes-Cummings model and describe its stationnary 
points, stable and unstable. We then explain and compute the classical monodromy when we loop 
around the unstable point. We also introduce a reduced system that will be important in subsequent 
considerations. In section [3] we define the quantum model and explain that the appearence of a 
default in the joint spectrum of the two commuting quantities is directly related to the monodromy 
phenomenon in the classical theory. In section [3] we perform the above reduction directly on 
the quantum system and study its Bohr-Sommerfeld quantization. The reduction procedure favors 
some natural coordinates which however introduce some subtleties in the semi classical quantization 
: there is a subprincipal symbol and moreover the integral of this subprincipal symbol on certain 
trajectories may have unexpected jumps. We explain this phenomenon. The result is that the Bohr- 
Sommerfeld quantization works very well everywhere, except for energies around the critical one. 
In section [5] we therefore turn to the semiclassical analysis of the system around the unstable 
stationnary point. We derive the singular Bohr-Sommerfeld rules. Finally in section [6] we apply 
the previous results to the calculation of the time evolution of the molecule formation rate, starting 
from the unstable state. For this we need the decomposition of the unstable state on the eigenstates 
basis. We find that only a few states contribute and there is a drastic reduction of the dimensionality 
of the relevant Hilbert space. To compute these coefficients we remark that it is enough to solve the 
time dependent Schrodinger equation for small time. We perform this analysis and we show that 
we can extend it by gluing it to the time dependant WKB wave function. This gives new insights 
on the old result of Bonifacio and Preparata. 

2 Classical One-spin system. 

2.1 Stationary points and their stability 

Hence, we consider the following Hamiltonian 

H = 2es^ + ujbh + g {bs~ + bs+) (1) 

Here = s^,s^ = zt is^ are spins variables, and 6, 6 is a harmonic oscillator. The Poisson 
brackets read 

K,/} = -e,fe,s^ {b,b} = i (2) 

Changing the sign of g amounts to changing (6, b) (—6, —b) which is a symplectic transformation. 
Rescaling the time we can assume g = 1- 
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The Poisson bracket of the spin variables is degenerate. To obtain a symplectic manifold, we 
fix the value of the Casimir function 

so that the corresponding phase space becomes the product of a sphere by a plane and has a total 
dimension 4. Let us write 

H = Ho + LoHi 

with 

Hi = bb + s'', Ho = 2ks'' + bs- +bs+, K = e-uj/2 (3) 

Clearly we have 

{Ho,Hi} = 

so that the system in integrabl^. Hence we can solve simultaneously the evolution equations 

dtJ = {Ho,f}, dtJ = {Hi,f} 

Once f{tQ,ti) is known, the solution of the equation of motion dtf = {H,f} = (dtQ +ujdt^)f is 
given by f{t,ujt). The equations of motion read 



dt^b = —ib 


dtob = 


—is~ 




dtb = 


—is — 


iojb 


(4) 


dty = o 


dtos' = 


i(6s~ — 


bs+) 


dts' = 


i(6s^ — 


bs+) 


(5) 




dtos^ = 


2iKs'^ - 


- 2ibs^ 


dts+ = 


2ies+ - 


- 2ibs' 


(6) 


= —is~ 


dtoS" = 


—2iKS~ 


+ 2i6s^ 


dts~ = 


—2ies~ 


+ 2^65"^ 


(7) 



The ti evolution is simply a simultaneous rotation around the z axis of the spin and a rotation of 
the same angle of the harmonic oscillator part: 

6(to,ti) = e-**i6(to), s±(to,ti) =e±'*^s±(to), s^(to, ti) = s"(to) (8) 

In physical terms, the conservation of Hi corresponds to the conservation of the total number 
of particles in the system, when this model (with N quantum spins 1/2) is used to describe the 
coherent dynamics between a Fermi gas and a condensate of molecules [13j . 
We wrote the full equations of motion to emphasize the fact that the points 

s± = 0, s'' = ±Scu b = b^ = Q (9) 

are all the critical points (or stationary points i.e. all time derivatives equal zero) of both time 
evolutions and ti and hence are very special. Any Hamiltonian, function of Hq and Hi, will have 
these two points among its critical points. However, it may have more. 

^This remains true for the Ai'-spin system. In the Integrable Community this model is known to be a hmiting case 
of the Gaudin model [14] 
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For instance when H = Hq + ujHi an additional family of stationnary points exist when cj^e^ < 

ci- 



si,. They are given by 
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their energy is given by 



S' = -1(4+^V) 



UJ 

The energies of the configurations Eq. ([9]) are 

E = ±2esci 

so that 

E-E' = -{sci±eujf > 

and we see that E' represents the degenerate ground states of H, breaking rotational invariance 
around the z axis. 

We are mostly interested in the configurations Eq. Q. To fix ideas we assume e < so that 
among these two configurations, the one with minimal energy is the spin up. It looks like being a 
minimum however it becomes unstable for some values of the parameters. To see it, we perform 
the analysis of the small fluctuations around this configuration. We assume that b, b, are first 
order and = Scie + 5s^ where e = ±1 according to whether the spin is up or down. Then 6s^ is 
determined by saying that the spin has length Sd 

Js"" = —s~s~^ 

This is of second order and is compatible with Eq. ([5]). The linearized equations of motion (with 
respect to H) are 

b = —iujb — is~ (11) 
s~ = -2ies~ + 2iscieb (12) 

and their complex conjugate. We look for eigenmodes in the form 

b{t) = b{0)e-^'^\ s- = s-(0)e~2^^* 

We get from Eq. (fT2]) 

Inserting into Eq. pi|) . we obtain the self-consistency equation for E: 

E = ^-'-^^ (13) 
2 2 E-e ^ ' 
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The discriminant of this second degree equation for E is — 2scie. If e = —1 it is positive and we 
have a local energy maximum. However if e = +1, the roots become complex when k? < 2s d- In 
that case one of the modes exponentially increases with time, i.e. the point is unstable. Since we 
are in a situation where k < 0, the transition occurs when k = —\J2sci- 




Figure 1: The image in phase space of the level set of the critical unstable point = 2kSci, Hi = Sd 
is a pinched two dimensional torus. 

The level set of the unstable point, i.e. the set of points in the four-dimensional phase space 
of the spin-oscillator system, which have the same values of Hq and Hi as the critical point Hq = 
2kSci,Hi = Sci has the topology of a pinched two-dimensional torus see Fig.[l]. This type of 
stationary point is known in the mathematical litterature as a focus- focus singularity \25\ [26l [27] . 
The above perturbation analysis shows that in the immediate vicinity of the critical point, the 
pinched torus has the shape of two cones that meet precisely at the critical point. One of these cones 
is associated to the unstable small perturbations, namely those which are exponentially amplified, 
whereas the other cone corresponds to perturbations which are exponentially attenuated. Of course, 
these two cones are connected, so that any initial condition located on the unstable cone gives rise 
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to a trajectory which reaches eventually the stable cone in a finite time. Note that this longitudinal 
motion from one cone to the other is superimposed to a spiraling motion around the closed cycle 
of the pinched torus which is generated by the action of Hi. 




Figure 2: The image in of the phase space by the map {Hi = m — Sci,Hq). The green dots 
are the critical points. The stable one is the point on the vertical axis (m = 0, —2kSci) and the 
unstable one is located at (m = 2sci,2kSci)- The upper and lower green boundaries are obtained 
as explained in the text. 



An important object is the image of phase space into under the map {Hi,Hq). It is shown 
in Fig. [2]. It is a convex domain in R^. The upper and lower green boundaries are obtained as 
follows. We set x = bb and 

Hi = X + = m — Scl, Max(0, m — 2sci) < x < m (14) 

where the parameter m is introduced to ease the comparison with the quantum case. Then we set 

b = ^e'^,b = ^e-'\ = - [s^Ye^'^ = y^{m - x){2sd - m + x)e'^'^ 

so that 



Hq = 2K{m — Scl — x) + 2\J x{m — x){2sci — m + x) cos{6 — ip) 



It follows that 



2K{m — Scl — x) — 2\J x{m — x){2sci — m + x) < Hq < 2K{m — Sd — x) + 2\J x{m — x){2sci — m + x) 
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To find the green curves in Fig. [2] we have to minimize the left hand side and maximize the right 
hand side in the above inequahties when x varies in the bounds given in Eq. (jl4p . 

In the unstable case, because an initial condition close to the critical point leads to a large 
trajectory along the pinched torus, it is interesting to know the complete solutions of the equations 
of motion. They are easily obtained as follows. Remark that 

(s^)2 = -(Is- -bs+f = -{bs- +bs+f + Abbs+s- 
= -(Ho- 2Ks^f + A{H, - s^){sl - {s^f) 

or 

(s^f = A{s^f - A{H, + + Ai^Ho - sl)s^ + (44^1 - Hi) 

Hence, satisfies an equation of the Weierstrass type and is solved by elliptic functions. Indeed, 
setting 

= \{Hi + k'^)-X, s' = iY 

the equation becomes 

y2 = _ - 53 (15) 

with 

52 = -{Hi - 3HoK + 2HiK^ + K^ + 3sli) 

53 = J-(_27iJ^ - 8^1^ + SQHqHik - 2^1^^ + 36i?oK^ - 24Fi«:^ - + nHisli - 36k^s^;) 
Therefore, the general solution of the one-spin system reads 

= \{Hi + K^)-p{it + a) (16) 

and 

bb = Hi-s^ = -{2Hi - K^) + p{it + a) (17) 
3 

where a is an integration constant and p{9) is the Weierstrass function associated to the curve 
Eq. (jlSp . Initial conditions can be chosen such that when t = 0, we start from a point intersecting 
the real axis Y = 0. This happens when a is half a period: 

a = or UJ2 or lj^ = —uji — uj2 

This general solution however is not very useful because the physics we are interested in lies 
on Liouville tori specified by very particular values of the conserved quantities Hq,Hi. For the 
configuration Eq. ([9]) with spin up, the values of the Hamiltonians are 

Ho = 2ksci, Hi = Sd 
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In that case we have 



and 



where we have defined 



1 \ ^ / 9 



= 2sci — > in the unstable case. (18) 



Hence, we are precisely in the case where the elliptic curve degenerates. Then the solution of the 
equations of motion is expressed in terms of trignometric functions 

X{t) = - n'^tanh^{n{t - to)) 

and 

-'"-^^"> = coshes -fa) 

This solution represents a single solitonic pulse centered at time to- When the initial condition 
is close but not identical to the unstable configuration with spin up and 6 = 6 = 0, this unique 
soliton is replaced by a periodic sequence of pulses that can be described in terms of the Weiertrass 
function as shown above. Note that in the context of cold fermionic atoms, bb represents the number 
of molecular bound-states that have been formed in the system. 

2.2 Normal form and monodromy 

The dynamics in the vicinity of the unstable equilibrium can also be visualized by an appropriate 
choice of canonical variables in which the quadratic parts in the expansions of Ho and Hi take a 
simple form. Let us introduce the angle u g]0, 7r/2[ such that |«;| = ^/2sci cos u and $7 = \/2i^ sin u. 
As we have shown, instability occurs when < 2sci, which garanties that u is real. Let us then 
define two complex coordinates Ag and Au by: 



A, \ 1 f e-^''/^ 



The classical Poisson brackets for these variables are: 



(20) 



{As, As} = {A„,i4 = (21) 
{As,Au} = {As,Au} = l (22) 

Here we have approximated the exact relation {s+,s~} = 2is^ by {s"^,s^} = 2iscu which is 
appropriate to capture the linearized flow near the unstable fixed point. After subtracting their 
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values at the critical point which do not influence the dynamics, the corresponding quadratic 
Hamiltonian Hq and global rotation generator Hi read: 



Hi 



2n^{AsAu) + kHi 



(23) 
(24) 



With the above Poisson brackets, the linearized equations of motion take the form: 



A 



■u 



■s 



(-17 - iK)A_ 
{Q, — iK)Au 



(25) 
(26) 



With these new variables, the pinched torus appears, in the neighborhoood of the unstable point, as 
the union of two planes intersecting transversally. These are defined by A^ = which corresponds 
to the stable branch, and by = which gives the unstable branch. The global rotations generated 
by Hi multiply both Ag and A^ by the same phase factor e~**^ . We may then visualize these two 
planes as two cones whose common tip is the critical point, as depicted on Fig. [H The stable cone 
is obtained from the half line where Ag is real and positive and Au = 0, after the action of all 
possible global rotations. A similar description holds for the unstable cone. As shown by Eq. (I19p . 
any trajectory starting form the unstable branch eventually reaches the stable one. This important 
property will play a crucial role below in the computation of the monodromy. 

For latter applications, specially in the quantum case, it is useful to write down explicitely Ag 
and Au in terms of their real and imaginary parts: 



This reproduces the above Poisson brackets, provided we set {Xi,Xj} = 0, {Pi,Pj} = 0, {Pi,Xj} = 
6ij for i,j £ {1,2}. The quadratic Hamitonian now reads: 



which is the standard normal form for the focus-focus singularity [25 ^ 1261 [27]. 

We are now in a position to define and compute the monodromy attached to a closed path 
around the unstable point in the {Hq,Hi) plane as in Fig. [2j Let us consider a one parameter 
family of regular invariant tori which are close to the pinched torus. This family can be described 
by a curve in the {Hq,Hi) plane, or equivalently, in the complex plane of the AgA^ variable. Let 
us specialize to the case of a closed loop. Its winding number around the critical value in the 
{Hq, Hi) plane is the same as the winding number of AgA^ around the origin in the complex plane. 
Let us consider now the path x ^ [Oj 27r] AgAu = rfe^^ where > is assumed to be small. 
To construct the monodromy, we need to define, for each value of x, two vector fields on the 




(27) 
(28) 



Ho 
Hi 



n{XiPi + X2P2) + kHi 
X1P2 - X2P1 



(29) 
(30) 
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corresponding torus, which are generated by x-dependent hnear combinations of Hq and Hi and 
whose flows are 27r-periodic. Furthermore we require that the periodic orbits of these two flows 
generate a cycle basis (61,62) in the two-dimensional homology of the torus. This basis evolves 
continuously as x increases from to 27r. The monodromy expresses the fact that {ei, 62) can turn 
into a different basis after one closed loop in the {Hq, Hi) plane. The discussion to follow has been 
to a large extent inspired by Cushman and Duistermat |27j . 

One of these two flows can be chosen as the 27r global rotation generated by Hi. The cor- 
responding orbits are circles which provide the first basic cycle 61. For the other one, we first 
emphasize that the fiow generated by Hq is not periodic in general. So, to get the other ba- 
sic cycle 62, we have to consider the fiow generated by an appropriate linear combination of Hq 
and Hi. To construct it, it is convenient to consider an initial condition (^^(0), ^^(O)) such that 
As{0)Au{0) = rfe^^ and |^s(0)| <C |^u(0)|, that is close to the unstable manifold of the pinched 
torus. Let us pick a small enough C > so that the linearized equations of motion are still acurate 
when |^<j(0)p + |74tj(0)p < C^, and let us also assume that 77 <C C- Using the global rotation 

2 

invariance of the dynamics, we may choose ^^(O) = C, and As{Q) = ^e~^^. If we let the system 
evolve starting from this initial condition, \As\ will first decrease and \Au\ will increase, so the 
trajectory becomes closer to the unstable torus, along its unstable manifold. After a finite time 
ti, the trajectory reappears in the neighborhood of the unstable fixed point, but now, near the 
stable manifold. This behavior can be seen either as a consequence of Eq. (llOh for the trajectories 
on the pinched torus and extended to nearby trajectories by a continuity argument, or can be 
checked directly from the explicit solution of the classical dynamics Eqs. (I16|17jl in terms of the 
Weierstrass function. We may thus choose ti such that = A very important property of 

As{ti) is that it has a well defined limit when 77 goes to zero, because in this limit, one recovers the 
trajectory on the pinched torus such that As{0) = and Au{0) = C- As a result, when rj is small 
enough, the argument of As{ti) in polar coordinates weakly depends on x and the winding number 
of As{ti) when x goes from to 27r vanishes. We then let the system evolves until time ^2 using the 
linearized flow. The time t2 is chosen in order to recover |At((t2)| = C- Using Eq. (j26p and the fact 
that AsAu is conserved, this gives t2 — ti = (2/57) ln(C/?7). Furthermore, from As{ti)Au{ti) = rfe^^ 
and Eq. ([26|) . we deduce: 



In general, we have no reason to expect that Au{t2) should be equal to ^u(O) so, to get a periodic 
flow on the torus, the evolution generated by Hq during the time ^2 has to be followed by a global 
rotation of angle given by: 



The sign of (3 reflects the fact that the Hi flow applied during the time (3 multiplies Ag and A^ 
by 6~*^. Note that because the flows associated to Hq and Hi commute, the composition of the 




(31) 




(32) 
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i^o-flow during time t2 and of the Hi flow along an angle /3 can also be viewed as the flow generated 
by the linear combination Hq + {(3/t2)Hi during time ^2- We have established the periodicity of 

2 

this flow for the orbit starting from ^u(O) = C, and ^s(O) = ^e~*'^. But any other orbit on the 
invariant torus charaterized by AsAu = rfe^^ can be deduced from this one by a global rotation, 
which commutes with the flow of Hq + {[3/t2)Hi. This establishes the periodicity of this latter flow, 
and therefore allows us to construct the other basic cycle 62 (x); which depends smoothly on x- 

Because the winding number of As{ti) vanishes provided r/ is small enough, we have the crucial 
relation: 

/5(27r) = /3(0) + 27r (33) 

In other words, when we follow by continuity the periodic flow generated by Hq + {(3/t2)Hi during 
time t2, the final flow is deduced from the initial one by a global rotation of 2tt. To summarize, as 
X increases smoothly from to 2tt, ei is left unchanged, whereas 62 evolves into 62 + 61. In this 
basis, the monodromy matrix of the focus-focus singularity is then: 

M=(J ;) (34) 

2.3 Reduced system. 

As we have discussed before, the common level set of Hq and Hi which contains the unstable critical 
point is a pinched torus. Since we are mostly interested in the time evolution of the oscillator energy 
bb, and because this quantity is invariant under the Hamiltonian action of Hi, it is natural to reduce 
the dynamics to the orbits of Hi. For an initial condition on the pinched torus, this amounts to 
discard the spiraling motion around the torus, and to concentrate on the longitudinal motion from 
the stable cone to the unstable one. In other words, this reduction procedure reduces the pinched 
torus into a curve which has a cusp at the critical point. This is illustrated by the thick curve on 
Fig. El 

An other description of this trajectory is given by the thick curve on Fig. [Hthat will be discussed 
below. In practice, the first thing to do to implement this reduction is to fix the value of Hi to its 
critical value 

Hi = bb + s^ = Sci (35) 

This defines a submanifold Ti-s^i ™ phase space. Since Hq and Hi Poisson commute we may consider 
the reduced system obtained by performing a Hamiltonian reduction by the one parameter group 
generated by Hi . The reduced phase space is the quotient of TCs^, by the flow generated by Hi (the 
trivial ti phases in Eq. ([8])). It is of dimension 2. Explicitly, we set: 

X = bb> 0^b = V^6^^ b = ^e-'^ 

From Eq. (|35p and from the condition (s^)^ + s^s^ = s^; we deduce 

= sd. - X, and = ^x{2sci - x)e^''^ (36) 



13 



Figure 3: The sphere represents the reduced phase-space obtained by fixing the conserved quantity 
Hi to its value Sd at the critical point. Here, this point corresponds to the north pole on the sphere. 
The red curve is the image of the pinched torus after the reduction procedure. It is obtained by 
setting the other integral of motion Hq to its critical value Hq = 2kScI- For this, we have used the 
expression of Hq given in Eq. ()37p . The singular nature of the pinched torus is reflected by the 
cusp of the red curve at the north pole. 



the reduced Hamiltonian Hq reads 



Notice that x and 6 — ip are invariant by the Hi flow, and they can be taken as coordinates on the 
reduced phase space. These coordinates are canonically conjugate as we easily see by writing the 
symplectic form: 



In the following, when talking about this reduced system, we will simply set = 0. These coordi- 
nates are very convenient but one has to be aware that they are singular at x = and x = 2sci- The 
whole segment x = 0, < < 27r should be identified to one point and similarly for the segment 
X = 2sci, < 61 < 27r. 

The equations of motion of the reduced system read 



Hq = 2kScI + 2x\j2sci — X cos{9 — if) — 2kx 



(37) 



CO = -iSb ASb + i — - A 5s 



6x A 5{9 - if) 



X 



2x\/2sci — X sin 9 




From this we see that the critical points are given by 



cos 9 



K 



z + l = 0, z = e" 



x = 
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In our case k is negative, so that the solutions of the quadratic equations will be defined as 



z± = = -e^^", n"^ = 2sci - /^^ < z/ < 7r/2 

These points exist precisely in the unstable regime < 2s d- They correspond to only one point 
in phase space whose energy is given by 

Hq = 2kScI = (38) 

hence they correspond to the unstable point. The image of the pinched critical torus after the 
reduction procedure is shown on Fig. [4] as the thick line connecting z+ and z- . Note that this line 
intersects the x = segment with a finite angle. This angle reflects precisely the pinching of the 
torus at the unstable critical point in the original four-dimensional phase space. 

Another critical point is obtained by setting 9 = tt. Then 

1 



x = 2sd-X'^, with X = -{-n + ^K? + &Sd) (39) 

Since < x < 2scu we should also have < < 2s d which is the case when k > —^/2sd^ that is 
to say in the unstable region. The energy is given by 

= 2KSd -2{X + k) {2sd - ^2) (40) 

We see from that equation that, in the unstable region, Hq — 2KSd is always negative and the 
configuration Eq. (j39l) represents the ground state. 



Finally, the energy maximum of this reduced system is obtained for 6 = and x = 2sd — 
with Y = (k + \//?"+~6)/3. On the phase portrait shown on Fig. HI one sees a curve which contains 
the point at x = 2sd and which looks like a separatrix. However, unlike the situation around x = 0, 
we notice that it does not correspond to any stationary point of the unreduced Hamiltonian. In 
fact it is simple to show that the curve is actually tangent to the vertical line x = 2s d- On the 
sphere, this is the trajectory passing through the south pole. 

3 Quantum One-spin system. 

3.1 Energy spectrum 

We now consider the quantum system. We set as in the classical case 

H = Ho + ioHi 
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Figure 4: The phase portrait corresponding to the Hamiltonian Hq, Eq. (I37p . in the {6, x) variables, 
which are natural coordinates for the reduced system. The thick red line corresponds to the critical 
separatrix Hq = 2kScU which is the image of the pinched critical torus after symplectic reduction. 
Note that the unstable critical point is mapped onto the vertical axis at x = 0, so that the critical 
separatrix is indeed a loop. The blue dots correspond to the minimal and maximal values of Hq 
when Hi = s^- 
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with 

Hq = {2e-uj)s' + bU- +hs'^ = 2ks' + bU- +hs'^, Hi = b'^b + s' 
where as before 

K = e-uj/2 (41) 

We impose the commutation relations 

[b, b^] = h, [s+, S-] = 2hs\ [s^ s^] = ±hs^ 
We assume that the spin acts on a spin-s representation. 

s^\m) = hm\m), s^|m) = h\/ s{s + 1) — m(m it 1) |m it 1), m = —s, —s + 1, • • • ,s — l,s 
where 2s is integer. Of course 

We still have 

[Ho,Hi]=0 

and the quantum system is integrable. The Hamiltonian H can be diagonalized by Bethe Ansatz, 
but we will not follow this path here. For related studies along this line for the Neumann model 
see |29j . Note that recently, the non equilibrium dynamics of a similar quantum integrable model 
(the Richardson model) has been studied by a combination of analytical and numerical tools. [30j 
Let: 

e„ = (26t)"|0) (s+)*^-"| - s), Sup(0, M - 2s) < n < M (42) 
For all these states one has: 

HiCn = h{M - s)en 

Since Hq commutes with Hi we can restrict Hq to the subspace spanned by the e^. Hence, we 
write: 

M 

^ ""WW 

n=Sup(0,M-2s) " 
where the norm IIct^N of tliG vector is given by: 

II, ||2_.2A/-no2n (2^)!(^-n)!^! 

" {2s-M + n)\ 

Using 

HoCn = 2hK{M - s- n)en + —{M - n){2s + I - M + n)e„+i + 2hnen-i 



17 



the Schrodinger equation 
becomes: 



ih'^ = n^^J(n + l)(2s + 1 - M + n)(M - n)p„+i + hh^Jn{2s - M + n)(M + 1 - n)p„_i 
+2hK{M — n — s)pn 

and the eigenvector equation Hq^ = reads: 

\/{n + l)(2s + 1 - M + n)(M - n)pn+i + y^n{2s - M + n){M + 1 - n)pn-i (43) 

+ 2hK{M -n- s)pn = Epn (44) 

In this basis, the Hamiltonian is represented by a symmetric Jacobi matrix and it is easy to 
diagonahze it numerically. Varying M = 0, 1, • • • we construct the lattice of the joint spectum of 
Hq,Hi. We see on Fig. [5] that the lattice has a defect located near the unstable classical point. 
This defect induces a quantum monodromy when we transport a cell of the lattice around it. The 
rule to transport a cell is as follows. We start with a cell and we choose the next cell in one of 
the four possible positions (east, west, north, south) in such a way that two edges of the new cell 
prolongate two sides of the original cell, and the two cells have a common edge. We apply these 
rules on a path which closes. If the path encloses the unstable point, the last cell is different from 
the initial cell. 

The precise form of this lattice defect can be related to the classical monodromy matrix M, as 
shown by Vu Ngoc [23]. Choosing locally a cycle basis Cj on the Arnold-Liouville tori is equivalent 
to specifying local action-angle coordinates ( Jj, (pi) such that is obtained from the periodic orbit 
generated by Jj. After one closed circuit in the base space of this fibration by tori, the basic cycles 
are changed into e[ = Y2j ^ji^j therefore, the new local action variables J- are deduced from 
the initial ones by: 

n 

J[ = Y,M,,J, (45) 
i=i 

Heuristically, the Bohr-Sommerfeld quantization principle may be viewed as the requirement that 
the quantum wave functions should be 2tt periodic in each of the n phase variables (pi, which 
implies that Jj should be an integer multiple of h. When h is small, the discrete set of common 
eigenvalues of the conserved operators has locally the appearance of a regular lattice, whose basis 
vectors vi,...,Vn are approximately given by: dJi{vj) = h6ij. Because of the monodromy, this lattice 
undergoes a smooth deformation as one moves around a critical value of the mutually commuting 
conserved quantities. After one complete turn, the basic lattice vectors are changed into 
which are also approximately given by: dJ[{v'-) = h6ij. Taking into account the definition of the 



ih 



'dt 
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Figure 5: The lattice of the joint spectrum of Hq and Hi (or equivalently the integer M) in the 
unstable regime. The green point represents the classical unstable point {Hi = Sd, or M = 2s, 
and Hq = 2kSci)- The parameters are h-^ = s = 10, and k = -0.5y/2^i. This fi gure demonstrates 
the non-trivial quantum monodromy in this system. For this purpose, we follow by continuity an 
elementary cell of this lattice along two closed paths. When the path encloses the image of the 
critical point in the (Hq, Hi) plane, we see that an initial square cell does not preserve its shape 
at the end of the cycle, but it evolves into a parallelogram with a tilt. In other words, the lattice 
of joint eigenvalues exhibits a dislocation located at the critical value of -f^o and Hi. 
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classical monodromy matrix in Eq (j45p . this leads to the semi-classical monodromy on the joint 
spectrum: 

n 

v\ = Y,{M-%v, (46) 
i=i 

The discussion above was performed in the action angle variables (Ji, J2)- In the {Hq,Hi) vari- 
ables the basis v,v' should be transformed by the Jacobian matrix of the mapping (Ji, J2) — > 

{Ho{Jl,J2),HiiJi,J2)). 

On Fig. m we see that after a clockwise turn around the singular value in the {Hi, Ho) plane 
(which corresponds to a positive winding in the {Hq,Hi) plane), the lattice vectors vm and vhq 
are transformed into vm — vhq and vhq respectively. This is in perfect agreement with the classical 
monodromy matrix computed in section 12.21 



When M = 2s, i.e. when Hi takes its critical value corresponding to the unstable point, 

HiCn = hs Cn 

the Schrodinger equation simplifies to: 

i^-^ = h2{n + l)-y/ (2s - n)pn+i + h^n^J (2s + 1 - n)pn-i + 2hK{s - n)pn (47) 

Setting 

X = hn, Sci = hs 

this equation becomes: 



^^^^^if^ = (a; + h)^/{2sci - x)p{x + h)+ x^/{2sci + h- x)p{x - H) + 2k{sci - x)p{x) (48) 



dp{x 

Introducing the shift operator 



e* p{x) = p{x + h) 

this can be rewritten as: 

dp{x) 



ih- 



dt 



y^{2sci - x)e'^x + xe"*^ V(2sd - x) + 2k{s^i - x) p{x) (49) 



and we recognize the quantum version of the classical reduced Hamiltonian Eq. (j37p . but with a 
very specific ordering of the operators. 

The eigenvector equation Hq^ = reads: 

3 , . n r . 3 



m{n + l)v^(2s-n)p„+i + mnyj {2s + I - n)pn^i + (2;ik(s - n) - E)pn = (50) 
The eigenvalues are shown in Figure O in the stable case and in the unstable case. 
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Figure 6: The eigenvalues of the reduced Hq in the stable regime on the left and in the unstable 
case on the right. The origin of the Energy axis is the classical critical energy Ec = 2hsn. In the 
unstable case the ground state energy lies below this line and close to the classical value given by 
Eq. (jlO]) . ( = s = 30, K = —2\J2sci (stable case), k = —Q.b\/2sci (unstable case)). 



4 Bohr-Sommerfeld quantization 

In this section we examine the standard semi-classical analysis of the one spin system. There are 
two complications as compared to the usual formulae. One comes from the fact that the phase 
space of the spin degrees of freedom is a 2-sphere and hence has a non trivial topology. The second 
one is related to the fact that the Schrodinger equation is a difference equation and its symbol has 
a subprincipal part. Let us recall the Schrodinger equation: 

{x + h)\/2sci — X p{x + h) + x\/2sci + h — X p{x — h) — 2kx p{x) = hep{x) (51) 

where 

he = E — 2hns, X = hn 

The first thing to do is to compute the Weyl symbol h(x, 9) of the Hamiltonian operator. In 
standard notations, it is defined by: 

{Hp){x) = ^J e~^n^--y^'h[^,9^p{y)d9dy 

It is straightforward to check that, for the Hamiltonian of eq. (l5ip 



h{x, 9) = 2yx+^ ^ 2sd + ^ - ^ cos 9 + 2k{sci - x) 

Expanding in h we find the principal and sub-principal symbols 

h{x,9) = h^{x,9) + hh\x,9) + --- 
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where: 

h^{x,9) = 2x\/2sci — X cos 9 + 2k{sci — x) 



hHx,e) = ( ^2^;;^^+ — =^= ) COS0 



We recognize, of course, that }{^{x^ 6) is the classical Hamiltonian of the system Eq. ()37p . In partic- 
ular the symplectic form reads: 

to = dx Ad9 

Note that the definition of the Weyl symbol we have used here is motivated by the simpler case 
where the classical phase space can be viewed as the cotangent bundle of a smooth manifold. In our 
case, the 6 coordinate is 27r periodic, and we are not strictly speaking dealing with the cotangent 
bundle over the line parametrized by the x coordinate. To justify this procedure, we may first 
imbed our system in a larger phase space, where both 9 and x run from — oo to cxd. Because the 
symbol h{x, 9) is 27r periodic in 6*, the unitary operator T associated to the 27r translation of 9 
commutes with the quantum Hamiltonian H associated to the symbol h{x,9). It is then possible 
to perform the semi-classical analysis in the enlarged Hilbert space and to project afterwards the 
states thus obtained on the subspace of 2tt periodic wave-functions, which are eigenvectors of T 
with the eigenvalue 1. Imposing the periodicity in 9 forces x to be an integer multiple of h. From 
Eq. (iST]) . we see that the physical subspace, spanned by state vectors \x = nh) with < n < 2s, is 
stable under the action of H. 

An alternative approach would be to use a quantization scheme, such as Berezin-Toeplitz quan- 
tization [31], which allows to work directly with a compact classical phase-space such as the sphere. 
When Hi takes its critical value corresponding to the unstable point, it is easy to show that the 
eigenvector equation (j50p can be cast into the Schrodinger equation for a pure spin s Hamiltonian 
given by: 

H^s = 2ks^ + s+Vsd - Sz + Vsd - Sz s~ (52) 

Starting from the known Toeplitz symbols of the basic spin operators, one may derive Bohr- 
Sommerfeld quantization rules directly from this Hamiltonian i^eff [32], but we shall not explore 
this further in this paper. 

Returning to our main discussion we now choose a 1-form 7 such that: 

j{X,,o) = -h\x,9) (53) 

where Xi^o is the Hamiltonian vector field associated to the Hamiltonian h^{x,9), which is tangent 
to the variety h^{x,9) = E. Using x as the coordinate on this manifold, we find: 



Xi^o = 2x\j2sci ~ X sin 9 dx 

Hence, we may choose: 

7 = - J ( - + 1^7^ r I cot9dx 

2 \x 2{2sd - a 
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Under these circumstances, the Bohr-Sommerfeld quantization conditions involve the form 7 and 
read (see e.g. j23|) 

where C{E) is the classical trajectory of energy E, a is the canonical 1-form oj = da, fJ.c{E) is 
the Maslov index of this trajectory and n is an integer. Note that the integral of 7 over C{E) is 
completely specified by the constraint Eg. ([53]) . 

We now come to the fact that the 2-sphere has a non trivial topology. A consequence is that 
the canonical 1-form a does not exist globally. In the coordinates x,6 we may choose 

a = X dO 

but this form is singular in the vicinity of the south pole where x = 2sci- Let us consider a closed 
path on the sphere parametrized by the segment at constant x, 9 running from to 2it. The integral 
of a around this path is 2ttx. When x goes to 2s d, the integral goes to 47rsci which is the total area 
of the sphere. But the limit path for x = 2sci, corresponds on the sphere to a trivial loop, fixed 
at the south pole. The consistency of Bohr-Sommerfeld quantization requires then iirsd to be an 
integer multiple of 2iTh. Hence we recover the quantization condition of the spin: 2s should be an 
integer. 

Before checking the validity of Eq. (l54p in our system, it is quite instructive to plot the integral 
of 7 along the curve C{E) as a function of E. This is shown on Fig. [7| Here, we see two singularities 
at the values of E where C{E) goes through either the north or the south pole. At the south 
pole Jq(^^^ 7 merely jumps by vr, whereas at the north pole, a jump of vr is superimposed onto a 
logarithmic divergence. We can easily find an asymptotic formula for this integral for energies close 
to the critical value Ec = 2sc/k: 

^ = S "^^t^' + 2(7r - 1.) - vr e{AE) + 0{AE, AE\og\AE\) (55) 

here AE = E - E^ and Q{AE) = 1 if AS > and if AS < 0. 

We wish now to show that these jumps are merely an effect of working with polar coordinates 
{x,9). It is instructive to consider first a quantum system with one degree of freedom, whose 
classical phase-space is the {q,p) plane. From any smooth function h{p,q), we define the quantum 
operator 0\Y,h by: 

Polar coordinates are introduced, at the classical level by: 

(p^q) = V2x{cos9,sm9) 
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Figure 7: The integral of the sub-dominant symbol as a function of energy (the red curve). We 
see the logarithmic divergence for the energy corresponding to the singular critical point, and the 
gap — vr at the energy corresponding to the trajectory passing through the south pole. The green 
curve is the asymptotic formula Eq. (j55p in which we have added the gap at the south pole to ease 
comparison and to show that it remains quite accurate even well beyond its expected range of 
validity. 
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This definition imphes that dp A dq = dx A d9 . At the quantum level, we start with the usual {p, q) 
operators. In the Hilbert space L'^{q), we introduce the eigenvector basis {\n)} of + Q^, such that 
{p^ + q'^)\n) = h{2n + l)\n). We may then set x = X^„>o To define 6, such that [9,x] = ih, 

it is natural to enlarge the Hilbert space, allowing n to take also negative integer values as well. 
Then we set exp(z^)|n) = |n + 1). From a\n) = ^/n\n — 1) and a^\n) = v^r)~+T|n + 1), we see 
that, in the physical subspace n > 0, we have a = h~^^'^ exp{—i9)x^^'^ and a+ = exp(i^). 
The Weyl symbols, in the (x, 9) variables, of these operators are h~^^'^{x + h/l)^^"^ exp(=pz^). Using 
p = {h/2y^'^{a + a+) and q = (/i/2)^/^i(a — a+), we deduce that the Weyl symbols of p and q in 
the {x,9) variables are respectively \/2x + h cos 9 and \/2x + h sm9. So the Weyl symbol is not 
invariant under the non-linear change of coordinates from {p,q) to {x,9). For linear functions of 
p and q the symbol in {x,9), denoted by hpoi{x,9) is obtained from h{p,q) by substituting p and 
q by their classical expressions as functions of {x,9) and then by shifting x into x + h/2. Such a 
simple rule does not hold for more complicated functions h{p,q). Nevertheless, one can show that 
the substitution of x ^ x + h/2 gives the symbol hpoi{x, 9) up to first order in h. Practically, this 
means that: 

hpo\{x,9) = h{\/2x + hcos9, \/2x + hsm9) mod Oih^) 

or, more explicitely: 

/ipoi(x, 9) = h{p, q) + ^ (P^^P^ ^) + '^^^P^ + 

where {p,q) = \/2x(cos ^, sin0). This has the following important consequence: even if h{p,q) has 
no subprincipal part (i.e. it does not depend on h), hpoi{x,9) does acquire a subprincipal part 
h^^i{x,9). Now, along any classical trajectory associated to h{p,q), we have: 

dh dh ■ 

p-Q^{p,q) + q-Q^{p,q) =pq-qp = ^x9, 

so h^^^{x, 6) = 9/2. This striking result implies that, when we work in polar coordinates, the 7 form 
associated to /ipoi(a;, 9) may be chosen as —d9/2. This immediately explains why Jq^^^^ 7 jumps by 
— TT when the /i(p, q) = E orbit crosses the origin of the polar coordinates, because then J(j(^^^ d9 
jumps from to 27r, see Fig. [51 Note that, in general, we expect h{p, q) to have also a subprincipal 
part h^{p,q). But this gives an additional term to h^oi(^'^) which is deduced from h^{p,q) simply 
by the classical change of variables from {p,q) to {x,9). This term has no reason to display any 
singularity when the classical trajectory goes through the origin. Remark that around the south 
pole, if we set 

X = 2sci — y, p = \/2ycos9, q = ^y2ysm9 

then we can expand 

ho = -2Sci{k + Scl) + K \ \p + ^^"^ ] + g2 
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Figure 8: On this figure we see that —\§d9 equals zero if the integral is taken over the small ellipse, 
while it is equal to — vr if it is taken over the big ellipse enclosing the origin. This explains why the 
integral of the subprincipal symbol in polar coordinates jumps by =b7r when the orbit crosses the 
origin of the polar coordinates. 
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and we see that our system is equivalent to a shifted harmonic oscihator. 

Note that on the sphere, polar coordinates have two singularities, at the north and at the south 
poles, which correspond to x = and x = 2sci respectively. To generalize the above analysis to the 
sphere, we may thus use two charts {p, q) and {p' , q') such that: 

{P)Q) = \/2x(cos 6*, sin 0) 
{p\q) = \/2^(cose',sin0') 
X + x' = 2sci 

e + 9' = 

We therefore get dp A dq = dx A d9 = dx' A dO' = dp' A dq' . It is interesting to note that the Weyl 
symbols hpoi{x,9) and /ipQi(x',^') are obtained from each other by the classical transformation 
from (x, 9) to {x' , 9'). On the other hand, h{p, q) and h'{p' , q') corresponding to the same quantum 
Hamiltonian in the physical subspace < x = 2sci — x' < 2sci do not coincide, because shifting 
x into X + h/2 amounts to shifting x' into x' — h/2 instead of x' + h/2. This implies that it is 
impossible to construct a quantum operator in the 2s + 1 dimensional Hilbert space of a spin s for 
which both subprincipal parts of h{p,q) and h'{p',q') would vanish. 

To formulate the Bohr-Sommerfeld rule, it is convenient to use the (p, q) coordinates when we 
consider a set of classical orbits which remain at a finite distance from the south pole. In these 
coordinates, the subleading term /(^(^^ 7 and the Maslov index are both continuous when C{E) 
crosses the north pole. Going back to (x, 9) coordinates, we see that the jump in /(^^^^^ 7 has to 
be compensated by a jump in the Maslov index. So the Bohr-Sommerfeld formula ()54p can be 
expressed in (x, 9) coordinates, provided these jumps in the Maslov index are taken into account. 

To check Eq. (l54|l we evaluate ^ Region) for the exact values and we define the defect 5n by: 

^ Region) = n + 5n 

Denoting 

I-i = TTT / "> ^0 = 7^ / 7 

27rft Jc{E) Jc{E) 

5n = I-i + Iq- {n + \), if Cn < 2?iKs 

6n = I-i + /o + ^ — + if 2/iKs < e„ < —2hKS 

5n= /-I + /o + 1 - (n + i), if en>-2hKS 

The excellent accuracy of this Bohr-Sommerfeld rule is clearly visible on Fig. [9l which shows 
the defect J„ as a function of n. We see that it remains small everywhere, and we emphasize that 
it does not exhibit any accident when e„ crosses the value —2kSci, corresponding to the south pole. 
On the other hand, something special happens around = 2kScI (the energy of the north pole), 
which cannot be attributed to the use of polar coordinates, but which reflects the crossing through 
the singular orbit associated to the pinched torus. A detailed description of the energy spectrum 
in the vicinity of the classical unstable point is the goal of next section. 



27 



Sn 
0.2 r 



0.1 - 



50 



100 



150 



200 



-0.1 - 



-0.2 L 



Figure 9: To test the quality of the Bohr-Sommerfeld quantization, we plot the default 5n as a 
function of n which labels the energy eigenstates. We see that the accuracy of the usual Bohr- 
Sommerfeld rule is very good away from the critical energy. In particular there is nothing special 
at the south pole. 

5 Generalized Bohr-Sommerfeld rule 

This region of the spectrum requires a different quantization rule [HI [T71 [18l [191 (SHI [Ml [M] because 
the classical motion near the turning point at Xmin is strongly affected by the presence of the 
unstable point at x = 0. Therefore, we first need to analyze the small x behavior of the energy 
eigenstates, when E is close to 2kScI- In the intermediate regime, between XminiE) and XmaxiE), a 
standard WKB analysis is quite acurate. Near XmaxiE), we have a regular turning point, which can 
be described, as usual, by an Airy function. Gluing the wave function of the eigenstates between 
these three different regimes will give us the generalized Bohr Sommerfeld quantization rules. 

5.1 Small X analysis 

Here, we shall recast the normal form around the singular point, studied in section [2?2] within the 
framework of the reduced system. For this, we assume n « 2s or x = hn « 2sci- In that 
approximation the Schrodinger equation Eq. (l50p becomes: 



Epn = hv2sh{n + l)p„4.i + hV2sh npn~i + 2/i/t(s — n)p,, 



'n 



(56) 



This equation is linear in n and can be solved by Laplace-Fourier transform. Letting: 



oo 




(57) 



n=0 
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the equation reads: 
with: 

Let us set: 



Hn 



(58) 



d 



d 



z = e 



de ^^dz 



the equation becomes: 



{V2^l{z + z ^) -2k)z—+V2^z 



5' 



Small 



Small 



We are looking for a solution which is analytic in a neighborhood of the origin in the z-plane. This 
can be written as: 



Small 



1 



z 



1 



where z± are the solutions of the second order equation already introduced before: 



z + 1 = 0, z± 



and 



A. 



1 e — K 1 „ , 



2 ' ' 2n 2 

Note that, when h is small, the leading term in f{E) is {E — Ec)/2hQ. It is quite instructive to 
express directly this wave-function in the 6 variable. Denoting 



x{9) = (cos 9 + cos u) "^/^ exp i 



/(^)log 



cos 



COS 



we have 



Small 
y Small 



(0) 



^r^man(^)^ i exp(-7r/(£;))x(^), 



— TT + l'<9<7r — v 
TT — V < 9 < IT 



Small 



(9) = -i exp{-Tr f{E))xi9), -tt<9<-'k + v 



Notice that this has precisely the expected form in the semi-classical limit h Q. Indeed, when E > 
Ec, the wave function is mostly confined to the interval —'k + v < 9 < -k — u, and it is exponentially 
small in the classically forbidden regions it — u < 9 < tt and — tt < 9 < —tt + u. Such exponentially 
small prefactor is reminiscent of the behavior of a tunneling amplitude. This is consistent with 
our treatment of 9 as an unbounded variable, subjected to a 27r-periodic Hamiltonian. Enforcing 



29 



integer values of n requires the wave-function $^™^"(^) to be 27r-periodic in 6. If were 
identically zero in the classically forbidden intervals, its Fourier transform would develop a tail for 
negative values of n. So these evanescent parts of the wave- function are required in order to ensure 
that the wave function belongs to the physical Hilbert space n > 0. Note that when E < Ec, the 
classically forbidden region becomes the intervals — tt < < — vr + and ir — < 9 < tt. 

In the classically allowed regions, the semi-classical wave function is expected to take the text- 
book form: 

^Smaii(^) = a(^)expi (^^p- + SiiO)^ (59) 

where the three functions a, Sq and Si take real values. These functions satisfy the standard 
equations [33], involving the principal and subprincipal symbols hP{x,9) = Ec + 2^/2sdx{cos9 + 
cosi^) and h^{x,9) = \^2sci cos 9: 

h^{S'oi9),9) = E (60) 

^{a\9)d,h\S',{9),9)) = (61) 

d,h^iS'o{9),9)S[{9) = -h\S'oi9),9) (62) 

From the above expressions for "^^^^^^(9), we can check that it has exactly this semi-classical form. 
Most likely, this occurs because $^™^"(0) is an eigenstate of a quantum Hamiltonian which is 
derived by symplectic reduction from the quadratic Hamiltonian of the normal form discussed in 
section 12.21 Quadratic Hamiltonians are particularly well-behaved with respect to the — > limit 
in the sense that the Bohr-Sommerfeld formula gives the exact energy spectrum, and also that the 
full quantum propagator obeys similar equations as ([50]) and (fUT|) . 

Let us now discuss the x = nh representation which is quite useful from the physical standpoint. 

^ ^ /co2^7^ ^ ' Jc,2i7: \ z+ ) \ z 

where Co is a small contour around the origin. When E is sufficiently far (in a sense to be precised 
below) from its critical value Ec-, we may use the saddle point approximation to evaluate \I'^™^"(x, e). 
This yields again the expected semi-classical form, namely: 

pSmaii(^) = 6(x)e±^"/^ exp -i + Wi{x)^ (63) 

where &(x), Wo{x) and Si{x) satisfy: 



/i°(x,Wo(x)) = E (64) 

dx 



^ {b\x)deh\x,W^{x))) = (65) 



deh^{x,WQ{x))Wi{x) = -h^{x,WQ{x)) (66) 
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Note that Wo{x) is the Legendre transform of So{9), that is Sq{6) = x, Wo(a;) = 6 and So{6) + 
Wo{x) = x0. It is interesting to write down exphcitely this semi classical wave-function when 
x 3> 1. It reads: 

pSf^aiw ) ^ B+{E)- 1 — + B(E)— J — 

where B_{E) = B+{E) and: 

B+{E) ,.,_,/2i_ Yo..n,M „ 4;iJ] E-E, 



B_{E) 



^.(-^/^)exp.(2/(^)log^ + ^) (67) 



We note that the phase factor diverges when E reaches E^. This behavior is induced by the 
contribution of the subprincipal symbol to the phase of the wave-function. The prefactor exp(i(i/ — 
7r/2)) in the above expression does not jump when E crosses the critical value Ec- This is consistent 
with the analysis of the previous section, because we do find a tt jump in the Maslov index that is 
compensated by a vr jump in the contribution of the subprincipal symbol. 

Let us now precise the validity domain of the stationary phase approximation. When we in- 
tegrate over 9 to compute the Fourier transform p^™'^^^(;r), we get an oscillating integral whose 
phase factor can be approximated by cxp(zS'q (0(x))(0 — 6'(x))^/(2^)) where 9{x) is one of the two 
saddle points defined implicitcly by Sq{9{x)) = x. This oscillating gaussian has a typical width 
{M"^) = h/\S'i{9{x))\. We find that: 



llsci (cos 9 -t- cos vY 



When E goes to E^ at fixed x, 9{x) goes to ib(7r — i^), because cos0 -|- cosz/ = E / {2y/2scix). So 
\Sq{9{x))\ ~ 2$7x^/|-E — Ec\. On the other hand, we have seen that the amplitude of ^'S°i^"(6l) 
diverges when 9 = ±(7r — i^). The distance 59 between 9{x) and the closest singularity goes like 
\E — Ec\/2^x. The stationary phase approximation holds as long as 9{x) is far enough from the 
singularities, that is if (A^^) << (69)'^, or equivalently, 

1^ - Ec\ » 2hn (68) 

When this condition is not fullfilled, the stationary phase approximation breaks down. The 
dominant contribution to p^™^"(x) comes from the vicinity of 9 = ±(71 — 1^). We have therefore to 
consider a different asymptotic regime, where H goes to while e = {E — Ec)/h remains fixed. In 
this case, the two exponents A+, A_ are fixed, and the only large parameter in the Fourier integral 
giving pSmaii^^;) is x/h. One then obtains: 

—i{n—u)x/h i{-K—v)x/h 

/--"(x) = A^{E) + A_iE)-—r-— = pf'^\x) +/_-^^x) (69) 

(x/n)2-*A^) {x/h)^+'f^^> 
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with A_{E) = A+{E) and: 

4±(|) ^ ^^iu-./2) - mE)) ^ ^^^^ ^ ^ 

This analysis first shows that, in spite of the breakdown of the stationary phase approximation, 
^Smaii^^^ is stiU given at large x by a sum of two semi-classical wave- functions associated to the same 
principal and subprincipal symbols and h^. This implies that there will be a perfect matching 
between and the WKB wave functions built at finite x from the full classical Hamiltonian. 

This will be discussed in more detail below. The only modification to the conventional WKB 
analysis lies in the phase factor A^/A- between the ongoing and outgoing amplitudes, which differs 
markedly from the semi-classical B^/B-. In particular, we see that the full quantum treatment at 
fixed e provides a regularization of the divergence coming from the subprincipal symbol. Indeed, 
the singular factor log{4hO,/\E — Ed) is replaced by the constant log (2 sin z/). Such a phenomenon 
has been demonstrated before for various models [201 128j . 

5.2 WKB analysis 

We return to the Schrodinger equation: 

{x + h)\/2sci — X p{x + h) + x\/2sd + h — X p{x — h) — 2kx p{x) = hep{x) (71) 

where: 

he = E -Ec 

We try to solve this equation by making the WKB Ansatz: 

Expanding in h, to order ^i*^ we find 

Xy/2Scl - X cos Wq- KX = (72) 

which is nothing but the Hamilton-Jacobi equation on the critical variety. It is identical to Eq. (|64p 
where h^{x, 9) is the complete principal symbol, and the energy is taken to be Ec- In this procedure, 
the energy difference he is viewed as a perturbation, which can be treated by adding the constant 
term — e to the complete subprincipal symbol h}{x,6). Alternatively, we may write the Hamilton- 
Jacobi equation as: 

e*^o=^(x), ^2(x)-^^^=z(x) + l = 

^J2Scl - X 

The solution of the quadratic equation is: 

z±{x) = {k ± i^/n-^-x) = ^ (73) 
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so that z±{x)\x=o = e^''^'^ Hence: 



(74) 



where the two signs refer to the two branches of the classical trajectory. Notice that they also 
correspond to the two determinations of the square root — X. 
At the next order in h the equation for b{x) is: 

d 



(b\x)xV^^-xj = (75) 

Therefore, we may choose b{x) = {xV^"^ — x)~^^'^. 

The correction Wi{x) to the phase function satisfies Eq. (j66p with h^{x, 9) replaced by h^{x, 9) — 
e, that is: 

deh\x,W^{x))W[{x) = -h\x,W(i{x)) + e 

W[{x) is the sum of two terms, W[{x) = —^{x) + e6P{x), where j{x)dx is simply the one-form 
7 defined by Eq. ([53]). evaluated on the critical trajectory, and 5f3{x) = {dghP{x,WQ{x))) ^. To 
understand better the meaning of 6(3{x), we start with the action integral Jq^^^-^ P, where /? = 9dx 
is closely related to the canonical 1-form a = xd9, because a + /? = d{x6). Note that, by contrast 
to Q, /3 is defined on the sphere only after choosing a determination of the longitude 9. Let us now 
study the variation of this action integral when E moves away from Ec and is changed into Ec + he. 
We have: 

P- f (5 = he I ^dx 



C(Ec+he) JCiEc) JC{Ec) dE 



But 9{x,E) satisfies h°{x,e{x,E)) = E, so: 



f)9 

deh\x,e{x,E)) — {x,E) = l 



and finally: 



^ - / p = he 5p{x)dx 

C{Ec+he) Jc{Ec) JC(E^) 



Explicitely, the equation for Wi{x) reads: 

W[{x) = ( ^ + (76) 

So we have: 

e ' = exp ±- log — ^=f=^ =F >^ log „ _ _ ^ (77) 
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We can now expand p^^^{x) when x is small. We find: 



Pr''i-)-^AY^%)-^^i^;^ (78) 



where: 



Notice that the phase factor is the sum of two terms which have simple geometrical interpretations: 



S+ = kQ + 2sciv = vrfi^ - f e{x)dx = [ a 

Jo Jc+ 



is the classical action computed on the upper half of the classical trajectory, and 

is the regularized integral of the subprincipal symbol j{x) = j{x) — e5fj{x) on the same trajectory: 

) ^{x)dx = lim / ^(x)dx + f(E) log — 

Jo ^A-'O yj^^ h J 

Similarly, setting a; = — |^| and expanding in jj- and we find to leading order: 

"""f^' - ^e*-'^l"'-{--«*)'^l"' (7.) 
The relevant wave function is of course a linear combination of p^^^{x). 
5.3 Airy function analysis 

Again, we start with the Schrodinger equation Eq. ([7T]) . We set: 

pAiry(^) ^ (-1)V'''''(2;), X = nh 

{x + h)y/2sd - X p^"^{x + h) + x^/2sd + fi-x p^"^{x -h)+ 2kxp^''^{x) = -hep^^'^ix) (80) 
and we expand x = fi^ + ^, keeping the terms linear in ^. Remembering that = —k, we find: 



de V 2av2 j dc 



/ d ^Airv 1 ft. ^ ^ A -Airv n 
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The unique solution which decreases exponentially in the classically forbiden region ^ > is pro- 
portional to the Airy function: 

^Airy(^) ^ e-^^AiTy{a^/^{^ - b)) 

where: 



b = h{ ^- -^'^^-^^ 1 , c=^(l 



When ^ < it behaves like: 

We expand in and One has: 
In that approximation e~'^^ ~ 1, and we get: 

5.4 Gluing the parts together 

The small x analysis gave: 



which is valid when x « 2sci- In an intermediate regime, both the small x and the WKB 
approximation are valid at the same time, as can be seen by comparing Eq. ()69p and Eq. (I78p . So 
we can glue these two wave- functions as follows: 

/(Small/' N /(Small/' \ 

We can now prolongate the pY^^{x) functions up to the region x ~ n'^ - \(,\. Recalling the 
asymptotic form of the WKB wave function in that region Eq. (I79p , we find 
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This has to be compatible with the Airy asymptotic formula which embodies the boundary condi- 
tions Eq. (jSip . Hence we find the condition: 



-1 



(82) 



This equation determines the energy parameter e and is the generalized Bohr-Sommerfeld condition. 
Given the fact that: 



Eq. (j82|) becomes: 



^ (2 ~ ^"20") 



rf) 



-2i[2^+(4s+l)i.+ 



2f^+2(2s+l)z.+ l^log 



Taking the logarithm, we find the quantization condition: 



^Smg(en) = (^'^ + ' " ^ ^) = ^KSd + hen 



(83) 



(84) 



where: 



^Sing{e) = -ilog 



^ (2 ^2^2) I 2'*^ 



+ 



h 



+ 2{2s + l)iy + 



K I 8f]3 
log 



h^/2Sd 



To test this condition, we can compute Sn = $(e^^^'^*) — 27r (n + ^) for the exact values of the 
energies. In Figure[Tni we plot 5„ as a function of n, using for <I>(e) both the usual Bohr-Sommerfeld 
function: 

1 



^ Regie 



h 



a + 



C{E) 



7 



C{E) 



and the function ^sing{e). Near the singularity, the singular Bohr-Sommerfeld condition is much 
more accurate. 

In Figure [TT] a typical example of the components of the eigenvectors is shown, comparing with 
the exact result obtained by direct diagonalization of the Jacobi matrix Eq. (jSOp and the various 
results corresponding to the different approximations: small x, WKB, and Airy. The agreement is 
very good. 

From Eq. (184p we can compute the level spacing between two successive energy levels. It is 
given by: 

- 27r _ Inn 

e-n+l — e-n — -, T — / X VOOj 



^'singi^n) log 
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Figure 10: To compare the usual and the singular Bohr-Sommerfeld rules, the default 5n is plotted 
as a function of n. The red dots are obtained by using the usual Bohr Sommerfeld rule, while 
the green dots are obtained by using $sing(e)- The dashed vertical lines represent the interval 
\E — Ec\ < 2hQ of validity of the singular Bohr-Sommerfeld rule. 




200 



Figure 11: The components of an eigenvector close to the critical level. The dots are the exact 
result. The red curve is the result of the small x analysis. The green curve is the WKB result and 
the blue curve is the Airy result. 
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where we defined 



i d ^ r(i 



(i + ix) 



This function has a sharp minimum located at x = where it takes the value — (7 + 2 log 2) where 
7 is Euler's constant, see Fig. [T2l Hence to leading order, the smallest energy spacing is: 



Figure 12: The energy level spacing e^+i — as a function the the state label k. This plot shows 
clearly the accumulation of energy levels in the vicinity of the critical point. 



Ae 



Inn 
log h\ 



A detailed analysis of the level spacing in the vicinity of the critical level has been given recently 
where it has been applied to the description of the long time dynamics of the system. 



6 Evolution of the oscillator energy 

Once the eigenvectors and eigenvalues are known, we can compute the time evolution of the oscil- 
lator energy: 

x{t) = (s|e*~6^6e~'~|s) 

Since the observable 6^6 commutes with Hi, its matrix elements between eigenstates of Hi with 
different eigenvalues vanish. Because |s) is an eigenvector of Hi with eigenvalue hs, we can therefore 
restrict ourselves to this eigenspace of Hi. The time evolution can be computed numerically by 
first decomposing the initial state on the eigenvector basis: 



(86) 
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so that the wave- function at time t is: 

\^{t)) = Y,Cne-'^-'/''\^{En)) 

n 

In the stable regime > 2s cu we get quite regular oscillations with a rather small amplitude as 
shown in Figure [T3l 

x(t) 



0.012 - 




Figure 13: The excitation number of the oscillator as function of time in the stable case. The green 
curve is the small time result Eq. (f98l) . {fr^ = s = 30, k = — 2\/2i^, Sd = 1). Notice that the 
vertical amplitude is about one hundred times smaller than in the unstable case (compare with 
Fig.dlD. 

By contrast, in the unstable regime < 2sci, it is energetically favorable to excite the oscillator, 
and we get a succession of well separated pulses shown in Figure [TH Note that the temporal 
succession of these pulses displays a rather well defined periodicity, but the fluctuations from one 
pulse to another are relatively large. The remaining part of this paper is an attempt to understand 
the main features of this evolution. 

On Figure [15] we show the coefficients \cn\- Only those eigenstates whose energy is close to the 
critical classical energy Ec contribute significantly. A good estimate of the energy width of the 

1 /2 

initial state |,s) is given by AE = {{s\H'^\s) - {s\H\s)^)^ . From Eq. §0^, we find: 

AE = hV2s^i = hQ/ sinv (87) 

Comparing with the criterion \AE\ ^ 2hil. for the validity of the stationary phase approximation, 
we see that most of the eigenstates which have a significant weight in the spectral decomposition 
of the initial state actually belong to the singular Bohr-Sommerfeld regime. 

An important consequence of this observation is that we can compute x{t) by considering only 
the few relevant states. The result is shown in Figure [16] for a spin s = 30. We retained only six 
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x(t) 

1.5 - , 




Figure 14: The excitation number of the oscihator as function of time in the unstable case. {h~^ = 
s = 30, n = —0.5^/2sci, Sci = !)• It clearly shows an aperidodic succession of pulses. The green 
curve is the semiclassical result Eq. (I102p . 



50 KM) 15() 200"' 

Figure 15: The coefficients of the expansion of the initial state \s) on the eigenvectors \^{En)). The 
red dots are the exact values. We see that only eigenvectors close to the critical level contribute 
significantly. The curve is given by Eq. (I89p . which is deduced from the knowledge of the short 
time evolution. The agreement between this approximate result and the exact one is remarkable, 
even in the immediate vicinity of the critical level. 
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states and superposed the result to the exact curve obtained by keeping the 61 states. We can 
hardly see any difference between the two curves. 




Figure 16: The excitation number of the oscillator as function of time for s = 30. The blue curve 
is the exact result, and the red curve is obtained by keeping only 6 energy eigenstates. 



We can compute the coefficients |c„| of the decomposition of the initial state \s). From Eq. (I86p . 
we have: 



e ft c„, r ^ 



e ft 



\c{E)\'^ p{E)dE 



where we have approximated the sum over discrete energy levels by an integral. Since we already 
know that the integral is concentrated around the critical energy Ec, the density p{E) can be 
computed from the singular Bohr-Sommerfeld phase eq.(| 



, ^, dn 1 d ^ , ^, Id 
= dE = 2^d^^^-^(^) = 2^hd-e 



On the other hand, since 



|eo| 



where cq is the state Eq. (j42p with the oscillator in its ground state, we get: 



(s|e * ft \s) 



2lT 



2lT 



^{e,t)de = po{t) 
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where ^{9,t) is the solution of the Schrodinger equation with initial condition ^{9,t)\t=o 
Inverting the Fourier transform in Eq. (|88p . we arrive at: 



\c{E)\'p{E) = — / e'-po{t)dt 



°° Et 



OO 



Taking for pQ{i) the small time expression Eq. (j96p . which is valid for times t < to — ^'~^r^-: 
expect to find an expression valid in the energy interval \E — Ec\ > 2^1 1 ^ i.e. for energies far 
enough from the critical energy. However because of the factor l/|log^| this interval covers most 
of the range Eq. ([871) where we expect c{E) to be substantially non zero. Notice also that | ^| is 
the order of magnitude of level spacing in the critical region, so that the interval where the use of 
the small time wave function is not legitimate contains at most a few levels. Explicitly, computing 
the Fourier transform of Eq. ()96p . we find: 

Icm'piE) = '~ Z,E, (89) 

where: 

a{E) = E - Ec - Hk 

The coefficients |c(-E')p computed with this formula are shown in Fig. [T5l The agreement with the 
exact coefficients is excellent. 

For longer time scales, we may also infer that the quantum dynamics exhibits the Bohr fre- 
quencies that can be deduced from the singular Bohr-Sommerfeld quantization rule (I84p . As shown 
before, the typical spacing between these energy levels is given by Ae ~ fiiglf' This corresponds 
to a fundamental period At ~ | log h\/^}. The fact that the energy levels are not equally spaced, as 
shown on Figure [12] is responsible for the aperiodic behavior on time scales larger than At. Sys- 
tematic procedures to analyze the long time behavior, starting from a nearly equidistant spectrum 
have been developed [35J, but we have not yet attempted to apply them to the present problem. 
Note that the quantum dynamics in the vicinity of classically unstable equilibria has received a lot 
of attention in recent years [Ml [37] . 

Instead, motivated by the analysis of the previous section, we now present a discussion of the 
time evolution in the semi-classical regime. As for the stationary levels, we show that we have to 
treat separately the initial evolution, when x is still small and the system is close to the unstable 
point, and the motion at later times, for which the usual WBK approach is quite reliable. As we 
demonstrate, this allows us to predict the appearance of pulses, together with their main period 
At. This generalizes the early study by Bonifacio and Preparata who focussed on the k = 
case. 
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6.1 Small time analysis 

Another expression for the mean oscihator energy is: 

2s 



i{t) = hY,n\Pn{t)\' 



n=0 

where Pnit) is the solution of Eq. (|47p with boundary condition: 

Pn{t)\t=0 = Sn,0 

For small time, Pnit) will be significantly different from zero only for small n. So in Eq. (j47p we 
may assume n « 2s. The equation becomes: 

— hv2sh{n + l)pn+i + hV2sh npn-1 + 2/ik(s — n)pn (90) 

As for the stationary case, this equation is linear in n and can be solved by Laplace-Fourier trans- 
form. We define: 

oo 

^{e,t) = J2Pn{t)e'''', ^{9,t)\t=o = l (91) 

n=0 

The time evolution is given by: 



= -{,/2^i{e"' + e^**^) - 2^)90^ - i(2/ts + V^ie (92) 

This equation can be solved by the method of characteristics. We introduce the function y(6) 
defined by: 

dy^ 1 

de V2sd(e'^ + e-'^) - 2k 

Explicitely: 

1 , e^^ + e^^ , ie ^ ( cosh fly 

yy^) = TTTT^Qg , -j,. ' or else e = ^= k + «J7 



217 e*^ e-^'" ^/2s^ V sinh Q.y 

The time dependent Schrodinger equation becomes: 

cosh. 

dt-^ + dy^ = -i(2s + 1)k^' + n — —r-^ 

smh Uy 



whose solution reads: 



^fy t) = -»«(2s+i)t sinhllj/ _ 

^ sinhl7(2/-t) '^^ 
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Imposing the initial condition Eq. (|9ip yields ^{y — t,0) = 1 and then: 



-m(2s+i)t sinhl^j/ 



sinh i}{y — t) 
Returning to the variable 6, we get: 

m{e, t) = e--(2^+i)* . „^ . ^^^^^ (93) 

As for the stationary case, we note that this wave-function has exactly the form dictated by 
semi-classical analysis, where the full symbol associated to Eqs. ([UU]) . is h{x,6) = 2y^2i^(x + 
^/2)(cos + cos v), after we have subtracted the energy Ec of the unstable point. The generalization 
of Eq. ([UD|) to the time dependent case is: 

h'^ideSoie, t),e) + dtSoie, t) = o (94) 

But since ^ is independent of at t = 0, 8080(6,1 = 0) = 0, which means classically that the 
trajectory begins at x = 0. But because of the form of h^, we get h^{dgSo{6,t = 0),0) = for any 
6, which implies that dtSo{6,t = 0) = 0. This shows that So{9,t) = identically, in agreement 
with the fact that a classical trajectory starting at x = stays there for ever! So the time evolution 
manifested in Eq. ()93p is of purely quantum nature, and it is encoded in the evolution of the 
amplitude a(9,t) and the subleading phase Si{9,t). Writing ^' = a(0, t) exp(f5'i(0, t)), the next 
order in h gives: 

c>,./i°(0, e)de^ + dt^ + ~ (a,/i°(0, 9)) 1^ ^ + ih\0, 9)^ = (95) 

A simple check shows this is exactly the same equation as (I92p without the term —i2KS^, due to 
the subtraction of the energy E^. Expanding in e*^ we find: 



^iu—Qt Q—iu+Qt 



(96) 



It is instructive to consider the large n limit. Then Pn{t) is proportional to exp[— ■|(«;(2sci + 
h)t + (tt — i^)x)]. This is the dominant phase factor for a WKB state of energy E = Ec which 
is concentrated on the outgoing branch of the critical classical trajectory. Note that there also 
appears a quantum correction to the energy, corresponding to a finite e = k. On the physical side, 
this is quite remarkable, because we have just seen that the short time evolution is driven by purely 
quantum fluctuations, formally described by the subprincipal symbol. Nevertheless, the subsequent 
evolution is quite close to the classical critical trajectory, because the energy distribution of the 
initial state is quite narrow, as we have discussed. 
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Further information is obtained by looking at the probabihty distribution of the number of 
emitted quantas: 

\Pn{t)\' - 



2sci sinh^ fit 



1^2 + 2sci s'mh^ nt 



02 + 2sci sinh^ Qt 
This is of the form: 

f)2 



2sc/sinh fit. 

Hence we have a thermal distribution with a time-dependent effective temperature. Such behavior 
is due to the strong entanglement between the spin and the oscillator. In fact we somehow expect 
that quantum effects, such as entanglement, are likely to be magnified in situations where there is 
an important qualitative difference between the classical and the quantum evolutions. 

It is now simple to compute the mean number of molecules produced in this small time regime. 
We find: 

, , ^ sinh^ rit ^ / 7T , , 

x{t) = 2hsci ^2 ' ^ = V2sd - ^2 (97) 

The small time approximation is valid as long as the number of molecules is small and will break 
after a time to such that x(to) — 1- To leading order in h, this time scale is: 

to --^ log/. 
Note that in the stable case, we change Q into iil and then: 

/ s ^ sin2 Q,t ^ rr, 

xit) = 2hsd—^, n = Vk^^2s^i (98) 

In that case we never leave the small n regime and this approximation remains valid even for large 
time as can be seen on Fig. [131 

6.2 Periodic Soliton Pulses 

When time increases, the approximation n « 1 is not valid anymore and we must take into 
account the exact quantum Hamiltonian. In the regime h « x = nh « 2sci, we can still perform 
a WKB approximation. The previous discussion has shown that the semi classical wave-function 
is concentrated on the classical orbit of with energy E = Ec + Hk, that is e = k. We set: 

p{x,t) = e-^(2s.,+;^)tg-i(H/o(x)WiW5(^^^) 

where 1^0(2;) is given by Eq. and Wi{x) by Eq. ([771) which e is replaced by it actual value 
K. The only source of time dependence, at this level of approximation, arises from the transport 
equation of the amplitude. Let us define the local velocity v{x) on the classical trajectory by: 



v{x) = -deh°{x,WQ{x)) = ±2xV^ 



X 
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Prom this field, we obtain the Hamiltonian fiow on the classical trajectory by solving the differential 

equation dx/dt = v{x). Let us denote by x{xo, t) the solution of this equation which starts from xq 
at t = 0. Likewise, we may reverse the flow and define a;o(x,t). It is convenient to introduce the 
function u{x) such that du/dx = —l/v{x). The solution of the fiow is then: 

u{xQ) = u{x)+t (99) 

For our problem, we have: 

^2 



i{x) = -j 



u{x) = — dx, or else x = 



2xV^'^ - X ' cosh^ nu 

The transport equation simply states the conservation of the local density b'^{x,t): 

dtib\x,t))+d^ib\x,t)vix))=0 

The evolution of this one dimensional conserved fluid is characterized by the invariance of b'^v along 
the flow: 

b'^{x,t)v{x) = b'^{xo{x,t),0)v{xo{x,t)) = B^{u{x)+t) (100) 

which simply results from the conservation law and the fact that the velocity fleld does not depend 
on time. 

Then we have: 

If we assume that Bo{u)^ is peaked around to, say Bq{u)'^ = d{u — to), we flnd: 

cosh^ Cl{t — to) ^ ^ 

which is just the classical expression. There is a simple way to evaluate the time to. In fact there 
must be a regime where the small time formula and this semiclassical formula should agree, i.e. 

h^sinh^ntc^ ^ , or n^e^^* ~ 4J^2e2"(*-*o) (103) 

cosh^ J^(t - to) 4J^2 ^ ^ 



and this gives the time to: 
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This reproduces well the position of the first peak. Its height is also well reproduced if instead of 
the crude approximation Bq{u)'^ = 5{u — to), we glue the wave functions given by the small time 
solution and the WKB solution at a gluing time tg-. 



B0{U + tg 



K 



sinh(r2n) 



■ exp 



cosh($7n)3 ^ \^ x{tg) cosh(il'u)2 
where the normalization constant K is given by the condition: 

f oo 




Bo(ufdu = 1 



and x{tg) is the value of x given by Eq. ()97p taken at time t = tg (see Fig. [T7|) . The height of the 



<n (t) > 




Figure 17: Time evolution of the mean oscillator energy x{t) for k = Sd = 1 and h = 1/300. 

The blue line corresponds to the exact solution, the red one to the semiclassical solution and the 
yellow line to the small time solution. 



first peak can be estimated from Eq. poip . The maximum is obtained when we have the largest 
overlap between 1/ cosh^ Qu and Bo{u + t)"^ . This is happens with a very good approximation, for 
t = to such that the maximum of Bq{u + to)^ is at zero: 

The previous equation allows us also to discuss the freedom in the choice of the gluing point 
tg. Indeed we expect that does not depend on tg as long as tg is restricted to a temporal domain 
where the short time solution and the semiclassical one overlap. If we plot to as a function of tg, 
we see that the curve we obtain develops a plateau, which means that the value of to we obtain 
chosing tg inside this plateau is independent of tg. 
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Figure 18: Time to of the first pulse as a function of the gluing time tg between the short time 
solution and the WKB solution. The existence of a plateau shows that gluing can be done in an 
unambiguous way. k = — 1/\/2, Sd = 1 and h = 10~^. We can recognize a plateau between tg = 1 
and tg = 4. 



Moreover, if we plot the value of to at the plateau, as a function of log(^), we recover with a very 
good accuracy the rough estimation of Eq. (jl04p . see Fig. [19l 

If we return to the de/ta-function approximation, it is clear that the situation reproduces itself 
after a time interval 2to- Hence we can write: 

y — ^- — — ^ n\n^[nit - to)\k^] 

where cn is a Jacobi elliptic function of modulus k'^ given approximately by: 

2hsci 



k^ ~ 1 



When K = 0, this is Bonifacio and Preparata's result 0. Clearly, it remains a challenge to account 
for the lack of exact periodicity in the pulse sequence, using the time-dependent approach. In this 
perspective, returning to the integrable structure may be interesting in order to obtain an effective 
model around the critical point in terms of particle-hole excitations of Bethe pseudo-particles. 



7 Conclusion 

Let us emphasize the main results of this work. First, the Jaynes-Cummings model exhibits, in a 
large region of its parameter space, an unstable fixed point which corresponds to the focus-focus 
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Figure 19: Time to of the first pulse as a function of h. The blue points are obtained by looking at 
the value of to at the plateau. The red line corresponds to the Eq. ()104p . k = — 1/\/2, Sd = 1- 



singularity of an integrable system with two degrees of freedom. As a result, there is a classical 
monodromy when one considers a loop which encircles the critical value in the {Hq,Hi) plane. At 
the quantum level, this phenomenon is manifested by a dislocation in the joint spectrum of Hq 
and Hi. We have then analyzed the eigen-subspace of Hi which corresponds to the critical point. 
The associated reduced phase-space is a sphere. We have shown how to perform a semi-classical 
analysis using the convenient but singular coordinates {x, 9) on this sphere. The main result here 
is that when the classical orbit crosses either the north or the south pole, where the longitude 6 is 
not well defined, the action integral associated to the subprincipal symbol jumps by itvr, and this is 
compensated by a simultaneous jump in the Maslov index. Most of the spectrum is well described 
by usual Bohr-Sommerfeld quantization rules, at the exception of typically |log?i| eigenstates in 
the vicinity of the critical energy, for which special Bohr-Sommerfeld rules have been obtained. 
Remarkably, the classical unstable equilibrium state, where the spin component is maximal and 
the oscillator is in its ground-state, has most of its weight on the subspace spanned by these singular 
semi-classical states. This fact explains rather well the three time scales observed in the evolution 
of the mean energy of the oscillator. At short time, this energy grows exponentially, reflecting 
the classical instability of the initial condition. At intermediate times, the energy of the oscillator 
exhibits a periodic sequence of pulses, which are well described by the classical motion along the 
pinched torus containing the unstable point. Finally, the delicate pattern of energy levels, which 
are not exactly equidistant, governs the aperiodic behavior observed for longer time scales. 

This work leaves many unsolved questions. One of them is to develop a detailed analytical 
description of the long time behavior. This should a priori be possible because, as we have seen, 
the initial state is a linear superposition of only a small number of energy eigenstates, for which 
the singular WKB analysis developed here provides an accurate modeling. Another interesting 
direction is the extension to several spins, which is physically relevant to the dynamics of cold 
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atom systems after a fast sweep through a Feshbach resonance. It would be interesting to discuss 
if the notion of monodromy can be generaUzed with more than two degrees of freedom. Finahy, 
it remains to see whether the qualitative features of the time evolution starting from the unstable 
state remain valid for an arbitrary number of spins, as may be conjectured from the structure of 
the quadratic normal form in the vicinity of the singularity. 
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